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Numerical simulations of the supernova (SN) neutrino self-induced flavor conversions, associated 
with the neutrino-neutrino interactions in the deepest stellar regions, have been typically carried 
out assuming the "bulb- model" . In this approximation, neutrinos are taken to be emitted half- 
isotropically by a common neutrinosphere. In the recent Ref. [l[ we have removed this assumption by 
introducing flavor-dependent angular distributions for SN neutrinos, as suggested by core-collapse 
simulations. We have found that in this case a novel multi-angle instability in the self-induced 
flavor transitions can arise. In this work we perform an extensive study of this effect, carrying out a 
linearized flavor stability analysis for different SN neutrino energy fluxes and angular distributions, 
in both normal and inverted neutrino mass hierarchy. We confirm that spectra of different v species 
which cross in angular space (where F Ve — F Vx and Fv c — Fv x ) present a significant enhancement of 
the flavor instability, and a shift of the onset of the flavor conversions at smaller radii with respect 
to the case of an isotropic neutrino emission. We also illustrate how a qualitative (and sometimes 
quantitative) understanding of the dynamics of these systems follows from a stability analysis. 

PACS numbers: 14.60. Pq, 97.60.Bw 



I. INTRODUCTION 

Flavor conversions of neutrinos emitted from core- 
collapse supernovae (SNe) represent a diagnostic tool to 
get crucial information about their mixing parameters 
and the SN dynamics in the deepest stellar regions 0. 
In particular, SNe are a unique laboratory to probe neu- 
trino oscillations in extreme conditions. Supernova neu- 
trinos can interact not only with the stellar medium via 
the Mikheyev-Smirnov-Wolfenstein (MSW) effect @, [|, 
but also with other neutrinos (v) and antineutrinos (v) 
as well. It was pointed out that large v densities in the 
deepest stellar regions can result in significant coherent 
v-v forward scatterings 0, @|. A few years ago it was 
discovered that v self-interactions can give rise to col- 
lective v flavor oscillations inside the SN 

EzHH (see 

for a recent review). The most important observational 
consequence of this collective behavior is a swap of the 
v e and v e spectra with the non-electron v x and v x spec- 
tra in certain energy ranges [ill EH ■ It has been argued 
that "spectral splits" at the edge of each swap interval 
would be observable in the high-statistics v signal from 
the next galactic SN, allowing to get crucial information 
about the unknown v mass ordering (see, e.g., [l3|). 

Self-induced oscillation effects crucially depend on the 
inner boundary conditions fixed for the further flavor evo- 
lution. Indeed, these non-linear flavor conversions are as- 
sociated to instabilities in the flavor space, that develop 
around the crossing points in the energy spectra of the 
different v species (where F„ c = F„ x and F Pc = ) [l^] ■ 
The number and the position of the crossing points de- 
pend on the ordering of the original SN v fluxes. Since 
this latter can change during the different post-bounce 
stages, significant temporal variations are expected in 
the pattern of the spectral splits [IH, [l4| HH- An im- 



portant layer of complication in this description is asso- 
ciated to the the current-current nature of the v-v weak- 
interaction Hamiltonian. This implies that the interac- 
tion energy between neutrinos of momenta p and q is 
proportional to (1 — v p • v q ) , where v p is the neutrino ve- 
locity 0, 0, [l(| . In a non- isotropic medium this velocity- 
dependent term would not average to zero, producing 
a different refractive index for neutrinos propagating on 
different trajectories. This is the origin of the so-called 
"multi-angle effects," which in some case can dramat- 
ically affect the development of the self-induced flavor 
conversions, producing a quick flavor decoherence [l7l.ll8j 
or suppressing flavor conversions otherwise possible for 
an a isotropic neutrino emission flil . [TBj . Moreover, it 
has been shown that the presence of the ordinary mat- 
ter background would cause a multi-angle suppression 
of the collective oscillations, when the matter density 
dominates over the neutrino one [l9j . This situation 
is expected to occur during the early times accretion 
phase [Hill- 



The characterization of multi-angle effects is then a 
key-ingredient to obtain a realistic description of the self- 
induced neutrino flavor conversions. In this context, it 
is expected that the v angular distributions at emission 
would play an important role in determining the v-v in- 
teraction strength. Until recently numerical simulations 
for the SN v flavor conversions have been based on the so- 
called "bulb model" (see, e.g, @,[Hl), where z/s of differ- 
ent species are considered as emitted "half-isotropically" 
(i.e. with all outward-moving angular modes equally oc- 
cupied and all the backward-moving modes empty) by a 
common spherical "neutrinosphere," in analogy with a 
blackbody emission. However, realistic supernova simu- 
lations show that v angular distributions at decoupling 
are far from being half-isotropic and, above all, are flavor- 
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dependent (see, e.g., [12, HH, HI]). The presence of non- 
trivial angular distributions was claimed in [13] to pro- 
duce a novel multi-angle instability in the self-induced 
flavor evolution of a toy model of v gas. However, that 
conclusion was based on an analysis performed with a 
small number of angular modes, making challenging to 
rely on this claim to infer conclusions for the realistic SN 
v case without a dedicated large-scale numerical study. 
Triggered also by this warning, in [l[ we performed nu- 
merical simulations of the self-induced flavor conversions 
with non-trival v angular distributions, finding remark- 
able effects on the flavor evolution. In particular, when 
flavor-dependent angular distributions lead to crossing 
points in the angular spectra of different v species, a 
new multi-angle instability can develop, in analogy to 
the known instability triggered by crossing points in the 
energy domain. We find cases in which this multi-angle 
instability can shift the onset of the flavor conversions 
toward low-radii and produce a smearing of the splitting 
features found with trivial v emission models. In order 
to achieve a semi-analytical understanding of the phe- 
nomenon, in that Letter we proposed to apply to our 
problem the linearized stability analysis recently worked 
out in 28]. By seeking an exponentially growing solu- 
tion in the eigenvalue equations associated with the lin- 
earized equations of motion for the neutrino ensemble, 
this method allowed us to determine the onset of the fla- 
vor conversions. In a specific scenario, we compared the 
growth of this mode in the half-isotropic case with differ- 
ent cases with non-trivial v angular distributions, finding 
a significant enhancement of the instability in these latter 
cases. 

The purpose of this follow-up work is to take a closer 
look to the multi-angle instability, triggered by non- 
trivial angular distributions. In particular we aim to use 
the stability analysis to perform an investigation on the 
dependence of this effect on the initial SN v flux ordering, 
and on the v mass hierarchy. Here is the plan of our work. 
In Sec. [TT] we introduce the setup for the flavor-stability 
analysis, describing the non-linear equations for the v 
flavor evolution in SNe, and the consistency equations 
coming from their linearization. In Sec. IIIII we present 
our models for the supernova neutrino emission for cases 
with a different flux ordering and angular distributions. 
For these different cases we show the results of the stabil- 
ity analysis in the two neutrino mass hierarchies. Finally 
in Sec. IIVI we comment on our results and we conclude. 



II. SETUP OF THE STABILITY ANALYSIS 

A. Equations of motion 

Our main goal is to perform a numerical stability anal- 
ysis for the cases studied in Therefore, we use the 
same setup followed in that paper. In particular, we work 
in a two-flavor oscillation scenario, associated to the at- 
mospheric mass-square difference Am 2 tm = 2 x 10 -3 cV 2 



and and with the small (matter suppressed) in-medium 
mixing = 10~ 3 . Since we aim at isolating the effect of 
the multi-angle instability, associated with v non-trivial 
angular distributions, we assume in the following that 
self-induced flavor conversions are not matter suppressed 
(as expected instead at t < 1 s after the core bounce), ig- 
noring the ordinary matter background in the equations 
of motion. Therefore, also when we will consider neu- 
trino spectra representative of the accretion phase (at 
t < 0.5 s after the core bounce) we will neglect the mat- 
ter term in order to isolate only the effect of the angular 
instability. However, it has been found in [22|, [23[ that 
taking realistic angular distributions during the accretion 
phase and a dominant matter term, self-induced flavor 
conversions would be completely inhibited. Conversely, 
we checked that the sub leading matter term during the 
cooling phase (at t>l s after the core bounce) would 
have a minor impact on the flavor instability associated 
with the non-trivial angular distributions. 

If we denote with i? r the angle of a given neutrino tra- 
jectory with respect to the radial direction, flux conserva- 
tion implies that that it will depend on the radial coordi- 
nate r even for straight line propagation (as we consider 
here). It is thus convenient to parameterize every an- 
gular mode in terms of its emission angle relative to 
the radial direction of the neutrinosphcre, that here we 
schematically fix at R = 10 km. For simplicity, we ne- 
glect possible residual scatterings that could affect z/s af- 
ter the neutrinosphcre, producing a small "neutrino halo" 
that would broaden the v angular distributions p9| . In- 
deed, it is expected that this effect may be relevant only 
at early times, where self-induced conversions would be 
matter suppressed [24|. For a half-isotropic distribution 
the occupation numbers are distributed as dn j d cos fin = 
const., or equivalently the radial fluxes are distributed as 
d<E> jd cos oc cos [l8| . A further simplification is ob- 
tained if one labels the different angular modes in terms 
of the variable u = sin 2 ^, as in [lj| HH- Note that 
for an half-isotropic emission at the neutrinosphere the 
v angular distribution of the radial fluxes is a box spec- 
trum in < u < 1, since dcos'&u/du oc (cosz?^)" 1 which 
cancels the cos fin dependence previously mentioned. At 
a radius r, the radial velocity of a mode with angular 
label u is v u>r = (1 - uR 2 /r 2 ) 1 / 2 []J]. Following [H, 
we write the equations of motion for the flux matrices 
<$>e,u a s function of the radial coordinate. The diagonal 
<&e,u elements are the ordinary number fluxes F l/a (E,u) 
integrated over a sphere of radius r. We normalize the 
flux matrices to the total V e number flux 7V Pe at the neu- 
trinosphcre. Conventionally, we use negative E and neg- 
ative number fluxes for anti-neutrinos. The off-diagonal 
elements, which are initially zero, carry a phase informa- 
tion due to flavor mixing. Then, the equations of motion 
read HG^l 



id r $ E ,u = [He.u^eA (1) 
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with the Hamiltonian [1, 0, H, [IJ 



where 



He,. — — 



1 M 2 



V2Gf 
Anr 2 



dE' / dv! 



1 - VuVu' 
V U V U > 



(2) 



in which we neglected matter effects. The matrix M 2 of 
neutrino mass-squares causes vacuum flavor oscillations 
and the term at second lines represents the v-v refrac- 
tive term. In particular, the factor proportional to the 
neutrino velocity v u>r in the v-v interaction term implies 
"multi-angle" effects for neutrinos moving on different 
trajectories 0-0, H(| ■ In order to properly simulate nu- 
merically this effect one needs to follow a large number 
[0(1O 3 )] of interacting v modes. 



B. Stability conditions 

In order to perform the stability analysis we closely fol- 
low the prescriptions presented in (28j and summarized 
in the following. Firstly, we switch to the frequency vari- 
able u) = Am 2 tm /2E so that E(oj) — \Am 2 tm /2oj\ and 
we introduce the neutrino flux difference distributions 
g u _ u = g(uj,u) defined as 



I Am? 



x 9H [F Vc {E(w),u) - F Vx (E(lj),u)] 



(3) 



2uj 2 



+0(-w) [F Vx (E(cj), u) - F Ve (E{w), «)] 



normalized to the total v e flux at the neutrinosphcre. 
Whenever referring to numerical values for the variables 
u),fi we shall implicitly quote them in km -1 , as appro- 
priate for the SN case. Note that g u>u is defined also for 
negative u, where it represents the difference of fluxes in 
the antineutrino sector in the opposite ordering. Then, 
we write the flux matrices in the form 1281 



Tr$„ 



9u 



(4) 



where Tr <f> UtU is conserved and then irrelevant for the fla- 
vor conversions, and the initial conditions for the "swap- 
ping matrix" in the second term on the right-hand side 
are s UfU = 1 and S U)U = 0. Self-induced flavor transi- 
tions start when the off-diagonal term S u u exponentially 
grows. 

In the small-amplitude limit |iS Wl «| <C 1, and at far 
distances from the neutrinosphere r>iJ, the linearized 
evolution equations for Suj,u m inverted mass hierarchy 
(IH, Aml tm < 0) assume the form [H| 

id r Su,u = (w + uefx)S UlU 

- y, du'duj'(u + u / )g LJ >,u>Su>, u > , (5) 



e = du dcu g^ 



(6) 



quantifies the "asymmetry" of the neutrino spectrum, 
normalized to the total v e number flux. The v-v interac- 
tion strength is given by 

V2G F N^ R 2 



4 7rr 2 2r 2 
3.5 x 10 5 / 



\10 52 erg/s 



15 McV 

(Eh.) 



R 



10 km 



One can write the solution of the linear differential equa- 
tion [Eq. J5|)] in the form S u , u = Qui,ue~ iar with complex 
frequency = 7 + in and eigenvector Q u u . A solution 
with k > would indicate an exponential increasing S^^, 
i.e. an instability. The solution of Eq. (J5)) can then be 
recast in the form of an eigenvalue equation for Q u ,u- 
Splitting this equation into its real and imaginary parts 
one arrives at two real equations that have to be satis- 
fied nn 



(Ji-m- 1 ) 2 

(Ji M 1 ) 



K( + J J 2 - K K 2 
J if 2 + K J 2 
2Ki 



(7) 



where 



J n = / dujdugu,uU n 



(u) - w ros ) 2 + K 2 



A" 



cLo dug u , u u r ' 



(lo - UJ rcs ) 2 + K 2 

and we introduced the resonant frequency 
w res (7) = 7 - uefi . 



(8) 



(9) 



A flavor instability is present whenever Eqs. ((7|) admit 
a solution (7, n). When an instability occurs, for a given 
angular mode Uq the function \Q UlU0 | is a Lorentzian [T^j, 
centered around w rcs , with a width n. Finally, we re- 
mind that the consistency equations in the normal mass 
hierarchy (NH) case (Am? tm > 0) are obtained simply 
changing the sign of g u . u in Eqs. ((HJ [28]. 



III. STABILITY ANALYSIS FOR DIFFERENT 
MODELS 

A. Models for supernova neutrino emission 

In order to perform our stability analysis we have to 
fix a neutrino emission model. Since our main goal is 
to explore in more detail the numerical findings pre- 
sented in [l[, in the following we refer to the cases dis- 
cussed in that work. We remark that energy and an- 
gular distributions of SN v's entering F Va (E, u) are not 
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FIG. 1: Spectrum g(oj) for the case A (upper panel) and C 
(lower panel). 
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FIG. 2: Difference of the energy-integrated angular spectra 
N Ve U Ve (it) — N Vm U Vx (u) for i/'s (left panels) and 77's (right 
panels) for the case A (upper panels) and C (lower panels). 
The different curves correspond to /3 e = f} x = 0.0 (dashed 
curves), f3 e — 1.0, f3 x = 1.5 (dotted curves), /3 e = 1.0, /3a; = 3.0 
(continuous curves), f3 e = 3.0, /3a, = 3.0 (dash-dotted curves). 



independent of each other. However, schematically we 
assume that the angular distributions are energy inde- 
pendent. Then, we can factorize the v flux of each flavor 
as F Va {E,u) = N Ua x ip Ua (E) x U Va (v). The v number 
N„ a = L Va I {E Va ) is expressed in terms of the v luminos- 
ity L„ Q and of the v average energy (E Ua ) of the different 
species. The function ip Va (E) is the normalized v energy 
spectrum (J dEip^ a (E) = 1) and U Va {v) is the normal- 
ized angular distribution (J duU Vct (u)=l). 

We parametrize the energy spectrum as in Ref. (3l| 



y(E) 



(1 + a) 



l+a 



E a 



T(l + a) 



cxp 



(l + a)E 



(10) 



where we fix the spectral parameter to a = 3 for all 
species. Following [l|, we fix the neutrino average ener- 
gies at 

«£„.), (E„ e ), (E Vm )) = (12, 15, 18) MeV . (11) 

Concerning the possible v flux ordering we consider two 
cases. As representative of the accretion phase (labelled 
as case A), we take 



N„ : : N„ = 1.50 : 1.00 : 0.62 



(12) 



corresponding to an asymmetry paramter e = 0.50 in 
Eq. (|6|). Instead, as representative of the cooling phase, 
dubbed case C, we choose 

N Ve : N Pa : N Vx = 1.13 : 1.00 : 1.33 , (13) 

giving e = 0.13. 



In Figure 1 we show the function — J du Quj,u 

for the 

case A (upper panel) and for the case C (lower panel). 
Both spectra present three crossing points, where = 0, 
so naively one would expect similar instability conditions 
in these two cases. However, as we will discuss in the 
following, the instabilities in these two cases are signifi- 
cantly different. 

Concerning the neutrino angular distributions, we use 
the simple toy model introduced in [l[ to capture the 
main deviations with respect to the half-isotropic bulb 
model, where XJ Va = 1 for all the v species. In particular, 
we choose forward-peaked distributions U„ a (u) cx (1 — 
k)/W 2 , For simplicity in the following we assume U u<! = 
Un and we refer to these cases: 



(A. 


= 0x = 0.0) 






= 1.0,/?.= 


1.5) 


G8e 


= i.o, Ac = 


3.0) 


We 


= 3.0 , A = 


3.0) 



In Figure [2] we we plot the difference N v JJ Va {u) — 
N Vx U Ujs (u) for i/s (left panels) and the analogous one 
for F's (right panels) for the case A (upper panels) and 
for the case C (lower panels) for the four (A>A) cases. 
The crossing points are given by the intersection of these 
functions with the horizontal lines at zero. In the half- 
isotropic case A = A = 0.0, the differences of angular 
spectra do not present any crossing point in the angular 
variable, while in all the other cases, the spectra present 
a crossing point at u — 1 where N v JJ Ve {u) = N^U^^u), 
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and all the angular distributions U Va vanish. In the 
case A the angular distributions do not present other 
crossing points at finite u. Also for the Vs in the 
/3 e = 1.0, f3 x = 3.0 case, the crossing point at finite u 
is so close to the edge of the spectrum at u = 0, that 
cannot be distinguished by it. Instead, in the case C the 
(anti)ncutrino angular spectra present also one crossing 
at finite < u < 1. In Qj, we associated the speed-up of 
the multi-angle instability with the presence or absence 
of crossing points in the energy-integrated angular spec- 
tra. In the following, we shall confirm this explanation 
by performing an extensive stability analysis the different 
cases. 



B. Case A 

The g u function for the case A (Fig. Q] upper panel) , 
presents three crossing points in the w variable. Accord- 
ing to what explained in [l2| and confirmed by the sta- 
bility analysis in (28|, the instabilities that trigger self- 
induced oscillations are associated to crossings with pos- 
itive slope in IH and to crossings with negative slope 
in NH. Therefore, one would have expected in this case 
a single self-induced spectral swap around lj = in 
normal hierarchy, and two spectral swaps in IH around 
the two crossings at finite cj. However, as observed 
in many numerical simulations for flux ordering with 
N v< , > N„ e > N Vx , the flavor evolution presents only 
a broad swap around ui = in IH and no conversion in 
NH. Indeed, as discussed in [l2j], a narrowly spaced triple 
crossing can superficially act like a single one at oj = 0. 

Hence, we present only the results for the IH case where 
the instability is expected to occur. Let us discuss first 
the half-isotropic case f3 e = f3 x = 0.0. In the case of 
a single-crossing spectrum, if the consistency equations 
[Eqs. ([7])] admit a solution, this is a unique (7, k) for a 
given fj,, or equivalently radial location in the SN [28j . 
In Figure [3] we show the (7, k) solution as a function of 
r (dashed curves for the half-isotropic case). The onset 
of the flavor conversions occurs when k starts to grow, 
i.e. at r ~ 90 km. In order to compare this prediction 
with the numerical results of the integration of the equa- 
tions of motion [Eqs. (JTJl-©], we represent in Figured] 
the integrated value of the z— component of the v polar- 
ization vector P z , that is related to the flavor content of 
the ensemble. The onset of the flavor conversions for the 
half-isotropic case (dashed curve) is in good agreement 
with expectations from the stability analysis. 

We pass now to analyze the cases with other neu- 
trino angular distributions presented in Sec. IIIA. Fig- 
ure [3] shows that different non-trivial cases do not lead 
to major changes in the behavior of the solution (7, n) of 
Eq. ([5]). In particular, no enhancement in the value of k 
is observable. This result is consistent with the numer- 
ical simulations of the flavor evolution (see Fig. 2]) that 
does not show any significant change in the non-trivial 
cases with respect to the half-isotropic one. The stabil- 
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FIG. 3: Inverted mass hierarchy, case A. Multi-angle eigen- 
values 7(r) and n(r). The different curves correspond to f3 e = 
j3 x = 0.0 (dashed curves), fj e = 1.0, f} a = 1.5 (dotted curves), 
/3 e = 1.0, f3 x = 3.0 (continuous curves), /3 e = 3.0, /3 X = 3.0 
(dash-dotted curves). 

ity analysis confirms that in absence of crossing points 
in the angular spectra, no new multi-angle instability is 
triggered by the non-trivial angular distributions alone. 
The only observable consequence of the non-trivial an- 
gular distributions is a shift in the onset in the flavor 
conversions toward smaller r. Since forward-peaked dis- 
tributions receive mostly contributions from small values 
of u, the resonance condition in Eq. <j9j> can be satisfied 
at larger /i, i.e. at smaller r with respect to the half- 
isotropic case. This is consistent with the expectation 
that the v-v strength is weaker for forward-peaked dis- 
tributions, making the system to become unstable earlier 
than in the half-isotropic case. In particular, we find as 
onset of the flavor conversions r — 87 km for the case 
with j3 e = 1.0, fix = 1-5, r = 88 km, for the case with 
(3 e = 1.0, (3 X = 3.0, and r = 78 km for j3 e = 3.0, (3 X = 3.0, 
consistently with what shown in Figure [U Finally, we 
also checked that no effect is triggered in NH by the pres- 
ence of non-trivial angular distributions. 

C. Case C 

The g u spectrum for the case C has been widely studied 
in the bulb model as representative of energy spectra 
with a flux ordering N Vx > N Ve > Np e and energy spectra 
with multiple crossing points leading to multiple splits 
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FIG. 4: Inverted mass hierarchy, case A. Radial evolution 
of the integrated ^-component of the neutrino polarization 
vector P z . The different curves correspond to fi e = fi x = 
0.0 (dashed curve), fi e = 1.0, fix = 1.5 (dotted curve), fi e = 
1.0, fix = 3.0 (continuous curve), fi e = 3.0, fi x = 3.0 (dash- 
dotted curve). 

around these points (see, e.g., [TH, Q, Il5j). Since we 
expect flavor instabilities in both the mass hierarchies, 
we discuss separately these two cases. 

1. Inverted mass hierarchy 

As discussed in [28j . with spectra presenting three 
crossing points the consistency equations [Eq. ([7])] ad- 
mit a pair of solutions (7+, k + ) and (7_, k_) (of course, 
if they admit a solution at all) for each value of [i, corre- 
sponding to negative and positive values of 7 respectively. 
We show in Fig. [5] and [6] the functions (7+(r), n+(r)) and 
(7_(r), K_(r)), respectively, for the different angular dis- 
tributions discussed before. A comparison between these 
two Figures reveals that ft_ is always smaller than k + , 
meaning that the instability associated with will be 
always sub-leading. Therefore, in the following it suffices 
to discuss only the (7+ , k+ ) solution. In the half- isotropic 
case /3 e = f3 x = (dashed curve) the n+ function presents 
a hump peaked around r ~ 100 km and connected with 
a long tail at smaller r around r ~ 75 km. This tail indi- 
cates that the system is in principle always unstable (also 
at very large values of fi) . However, at large y, and 7+ , the 
unstable frequency modes satisfying the resonance con- 
dition [Eq. (J9))] arc in the infrared region (i.e. |w| > 4), 
where the g u spectrum is strongly suppressed. Therefore, 
they have no impact for the flavor conversions. The onset 
of the flavor conversions is then given by the connection 
between the hump and the tail at r = 75 km. This re- 
sult is in agreement with what shown in the numerical 
calculation of P z , shown in Fig. (dashed curve) . 

When considering non-trivial angular distributions, 
the long tail in the k+ function observed in the half- 
isotropic case at small r now disappears. However, the 
presence of a non-isotropic angular distribution is not 
enough to produce an enhancement in the value of k+. 
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FIG. 5: Inverted mass hierarchy, case C. Multi-angle eigen- 
values 7+(r) and ft+(r). The different curves correspond 
to fi e = fi x = 0.0 (dashed curves), fi e = 1.0, Ac = 1-5 
(dotted curves), fi e = 1.0, fix ~ 3.0 (continuous curves), 
fi e = 3.0, fix = 3.0 (dash-dotted curves). 



Conversely, in the flavor blind case f3 e = (3 X = 3.0 
(dot-dashed curve) the instability is suppressed with re- 
spect to what has been seen in the half-isotropic case. 
Indeed, around the resonance energies [Eq. ([9|)], the 
integrals J„ and K n in Eq. ([7]) are proportional to 
j du dwg^.u u n k~ 2 . Since with forward-peaked distri- 
butions these integral receive mostly contributions from 
small u, in order to satisfy the consistency equations 
[Eq. (JT])] at a given fi, the corresponding k + has to be 
smaller than in the half-isotropic case. In this case, fla- 
vor conversions start at r = 80 km. 

A significant enhancement of the multi-angle instabil- 
ity occurs when v angular spectra exhibit crossing points 
in it, as can be seen in the case with (3 e = 1.0, (3 X = 1.5 
(dotted curve) and even more in the one with /3 e = 
1.0, f3 x = 3.0 (continuous curve). In both cases the peak 
in k+ increases with respect to the half-isotropic case 
and also the hump broadens toward smaller r. Since k + 
reaches a higher peak value, the width of the Lorentzian 
around an unstable frequency mode for a given angle u 
would be broad, implying a speed-up in the transitions. 
In particular, in the case of (3 e = 1-0, f3 x = 1.5 flavor con- 
versions start around r ~ 70 km, while in the case with 
p e = 1.0, j3 x — 3.0 around r ~ 55 km, in agreement with 
the numerical results shown in Fig. [7] 

In conclusion, we remark that SN v spectra with 
multiple crossing points in the energy domain can ex- 
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FIG. 6: Inverted mass hierarchy, case C. Multi-angle eigen- 
values 7-(r) and K-(r). The different curves correspond 
to /3 e = fi x = 0.0 (dashed curves), fi e = 1.0, fix — 1.5 
(dotted curves), fi e = 1.0, fix — 3.0 (continuous curves), 
fi e = 3.0, fix — 3.0 (dash-dotted curves). 



FIG. 8: Normal mass hierarchy, case C. Multi-angle eigen- 
values 7+(r) and K+(r). The different curves correspond 
to fi e = fix = 0.0 (dashed curves), fi e = 1.0, ft, = 1.5 
(dotted curves), fi e = 1.0, fix ~ 3.0 (continuous curves), 
fie = 3.0, fix = 3.0 (dash-dotted curves). 
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FIG. 7: Inverted mass hierarchy, case C. Radial evolution 
of the integrated ^-component of the neutrino polarization 
vector P z . The different curves correspond to fi e — fi x — 
0.0 (dashed curve), fi e = 1.0, fix = 1.5 (dotted curve), fi e = 
1.0, fix = 3.0 (continuous curve), fi e = 3.0, fix ~ 3.0 (dash- 
dotted curve). 



hibit flavor conversions at very large /i in a single- angle 
scheme [H]]. Multi-angle effects in a half- isotropic sit- 
uation s upp ress the instability at large neutrino densi- 
ties [l4l . |32| |. However, the stability analysis and our nu- 
merical simulations tell us that this stabilization effect 
in not a generic situation. Indeed, allowing for crossing 



points also in the angular spectra, flavor conversions can 
start much earlier than in the half-isotropic case. 



2. Normal mass hierarchy 

We finally pass to analyze the normal mass hierarchy 
case. Also here, we expect that the consistency equations 
[Eq. 0] admit a pair of solutions (7+, k+) and (7-, k_) 
for each value of /Lt, that we represent in function of r in 
Fig.[5]and[5]rcspcctivcly. The two imaginary parts of the 
eigenvalues, k + and k_ can assume comparable values, 
differently from the IH case. Whether k + or k_ plays the 
dominant role in triggering the flavor conversions will de- 
pend on the specific case. Starting with the half-isotropic 
case (fi e = fix = 0.0, dashed curves) we see that both 
the functions k + and k_ present long tails extending at 
small r. However, k+ is always smaller than 0.2, while 
k_ reaches a peak value ~ 0.6 at r ~ 105 km. Therefore, 
the onset of the flavor conversions will be associated to 
K-, and as usual is determined by the connection of the 
hump with the tail of the function, i.e. at r ~ 80 km, 
in agreement with the numerical evolution of P z shown 
in Fig. [10] (dashed curve for the fi e = fi x = 0.0 case). In 
the flavor-blind case with fi e = fi x = 3.0 (dash-dotted 
curves) the dominant instability is again the one asso- 
ciated with k_ that can reach as peak value ~ 0.5 at 
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FIG. 9: Normal mass hierarchy, case C. Multi-angle eigen- 
values 7-(r) and K-(r). The different curves correspond 
to /3 e = fi x = 0.0 (dashed curves), /3 e — 1.0, /3 X — 1.5 
(dotted curves), f3 e = 1.0, /3 X — 3.0 (continuous curves), 
/3 e = 3.0, f3 x — 3.0 (dash-dotted curves). 



slight larger radius, i.e. r ~ 90 km. This delay is prob- 
ably due to the time needed by the instability to grow 
significantly after it starts. This is related to the value 
of k_ that is relatively small in this case. In the case 
f3 e = 1.0, fi x = 1.5 (dotted curves) k_ is further reduced 
with respect to the previous cases, while n + is enhanced. 
Both the functions have a similar peak ~ 0.35, but «+ 
has a longer tail with a value close to the peak one, while 
k_ abruptly drops at zero for r ~ 95 km. Therefore, 
in this case the onset is determined by k + . Since the 
instable frequencies, determined by the resonance condi- 
tion [Eq. ©] would be in the infrared range (i.e. oj > 4) 
for fj, > 50, the onset of relevant flavor conversions cor- 
responds to this value of /x, i.e. r ~ 65 km, as visible 
in the numerical evolution (dotted curve). Finally, in 
the case /3 e = 1.0, fi x — 3.0 (continuous curve), the in- 
stability associated with k_ disappears, while the one 
associated with k+ is further enhanced, reaching as peak 
value ~ 0.4. This time the resonance condition reaches 
infrared frequency modes for [i > 27, corresponding to an 
onset of the flavor conversions at r < 77 km. Note that in 
this case the numerical evolution (continuous curve) the 
P z starts to vary earlier, around r = 55 km, but it rises 
significantly only at r > 70 km, as predicted by the sta- 
bility analysis. This can be interpreted as follows: in the 
latter two cases, since the K- function presents a rather 
flat behavior in /i, the determination of the onset of the 
flavor conversions in less accurate. Indeed, also the nu- 
merical evolution in Fig.flOlshows that the P z -component 
presents a slower rise time at the onset of the conversions. 
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FIG. 10: Normal mass hierarchy, case C. Radial evolution 
of the integrated ^-component of the neutrino polarization 
vector P z . The different curves correspond to /3 e — fix — 
0.0 (dashed curve), /3 e = 1.0, f3 m = 1.5 (dotted curve), /3 e = 
1.0, p x = 3.0 (continuous curve), /9 e = 3.0, fi x = 3.0 (dash- 
dotted curve). 



r ~ 105 km. We note that in this case K-(r) has no tail, 
terminating abruptly at r ~ 80 km, while n+(r) extends 
toward lower r reaching a almost constant value ~ 0.2. 
The predicted onset of the flavor conversions can then be 
located at r ~ 80 km. We see the corresponding numer- 
ical evolution in Fig. [TO] (dash-dotted curve) starts at a 



IV. CONCLUSIONS 

Self-induced flavor conversions for SN neutrinos are as- 
sociated to instabilities in the flavor space. Recently, a 
linearized stability analysis of the SN v equations of mo- 
tion has been proposed as diagnostic tool to track the 
emergence of these instabilities [28j . In our work we have 
applied this technique to study the multi-angle instabil- 
ity associated with non-trivial neutrino angular distribu- 
tions, extending the results presented in [l| to different 
SN v flux orderings and v mass hierarchies. We confirm 
that enhanced instabilities can occur, when the angular 
spectra present crossing points. Also, the onset of the fla- 
vor conversions typically shifts toward small radii, with 
respect to the half-isotropic case. 

We cannot but emphasize that neutrino angular distri- 
butions should be taken into account in supernova neu- 
trino phenomenological studies. For example, as pointed 
out in [l[, one can expect a smearing of the splitting 
features observed in the half-isotropic case, with result- 
ing v fluxes showing less significant spectral differences. 
This tendency toward spectral equalization would chal- 
lenge the detection of further oscillation signatures, like 
the ones associated with the Earth crossing of SN v's 
(see, e.g., [H, H3|)- On the other hand, a sufficiently 
low-radii onset on the flavor conversions might imply 
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an interesting impact on the r-process nucleosynthesis in 
SNe (35)] . Eventually, if v oscillations develop too close to 
the ncutrinosphere, they might invalidate the v transport 
paradigm in SNe that ignores v conversions. 

It is clear that future, state-of-the-art predictions of 
the neutrino flavor evolution in SNe would require as 
accurate as possible input from hydrodynamical simu- 
lations, providing the flavor, angular, energy and time 
structure of the fluxes. Here we showed that a prelim- 
inary flavor stability analysis presents some key advan- 
tages over a brute force numerical integration. While 
such a perturbative analysis cannot be used to deduce the 
eventual fate of the neutrino ensemble in flavor space, it 
can be useful to identify the situations where angular dis- 
tributions play a role in the flavor evolution. Once these 
cases are identified, it would be mandatory to perform 
large scale numerical simulations of the complete flavor 



evolution in order to determine the final oscillated neu- 
trino spectra. With current level of sophistication in SN 
simulations, we believe that such studies will be manda- 
tory to achieve a better characterization of the non-linear 
neutrino flavor evolution during a stellar gravitational 
collapse and a more realistic description of the neutrino 
signal from a future galactic SN. 
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